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1 Introduction 



Since Stoler, Saleh and Teich proposed the binomial states (BS) in 1985 |T[], so-called 
intermediate states which interpolate between some fundamental states such as number 
states, coherent and squeezed states and phase states have attracted much attention 0. 
The BS are finite linear combinations of number states 

M 
n=0 

where M is a non-negative integer, rj is a real probability (0 < r/ < 1) and \n) is a 
number state of the radiation field. The photon number distribution is clearly a binomial 
distribution, whence the name binomial state. The BS are intermediate number-coherent 
states in the sense that they reduce to number and coherent states in different limits 

r |M), r^->l, 

|r/,M)-^ |0), (1.2) 
[ \a), 7] ^ 0, M ^ oo, r]M = a'^. 

The BS also admit an eigenvalue definition 

(^y/^N + ^1 - vVM - Na^ \r], M) = ,/T]M\7], M), (1.3) 

where a, a'^ and are the annihilation, creation and the number operators, respectively. 
The algebra involved is the su(2) algebra (Holstein-Primakoff realization 0]) 

J+ = VM - Na, J- = aWM - N, = — - N, (1.4) 

2 

and in the present case the coherent state limit is essentially the contraction of su(2) to 
the Heisenberg-Weyl algebra generated by a^, a and 1. 

Since number and coherent states are eigenstates of the number operator and the 
annihilation operator a, respectively, it would seem more natural that, to define states 
interpolating between number and coherent states, we consider the eigenvalue equation 
of a linear combination of and a itself (not J"*"), namely, 

(^^N+^/l^a^ \\r],(3) = P\\r],(3). (1.5) 

Here < < 1 as before and (3 is an eigenvalue which will be determined not only from 
the eigenvalue equation Eq. (|1.5| ) but also by a physical requirement (see Sec. 2). 

In this paper we study the states ||?7, f3) and their various properties. We find that for 
f3 = y/rjM (M a non-negative integer), the solutions to Eq. (|1.5|) are indeed intermediate 



n 



M-n 



1.11 
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states which interpolate between number and coherent states. We also find that these 
states are closely related to the photon-added coherent states proposed by Agarwal 
and Tara The properties of this new state, such as their sub-Poissonian statistics, 
antibunching effects and squeezing effects, as well as their quasi-probability distributions 
(the Q and Wigner functions), are studied in detail. Although coherent and number 
states are not squeezed, the new interpolating states are squeezed, and exhibit highly 
nonclassical behavior. We also propose a scheme to produce these intermediate states 
in a cavity. 

The intermediate number-coherent states are of particular interest in their interac- 
tion with atomic systems. In the context of the Jaynes-Cummings (JC) model, atomic 
population inversion exhibits two completely different phenomena: Rabi oscillation when 
the field is initially prepared in a number state; and periodic collapse and revival when 
the field is initially prepared in a coherent state. We naturally expect that the states 
proposed in this paper will present phenomena intermediate between Rabi oscillation 
and periodic collapse-revival, given that the initial state of the field is in an intermedi- 
ate state. In Sec. 6 we study the interaction of the states with the atomic systems based 
on the two-photon JC model and we indeed observe this intermediate behavior. We 
also give an analytic derivation of the approximate photon number distribution of the 
field and find it exhibits strong oscillation at r = 7r/4,37r/4 (r the scaled time). This 
phenomena is explained physically in terms of the entropy and Q-function of the field. 

2 New intermediate number-coherent states 

In this section we solve the eigenvalue equation Eq.( |1.5| ), discuss the relation of the states 
(1.5) to photon-added coherent states and study the limit to number and coherent states. 

2.1 Solutions 

Expanding the state \\rj,l3) in number states 

oo 

h,/3) = ^C„|n) (2.1) 

n=0 

inserting it into Eq. ( |1 . 5|) and comparing the two sides of the equation, we find the 
solution of the eigenvalue equation Eq. (pr5|) 

h,l3)=Col^ ^ |n), (2.2) 

n=o (VI -^) vn! 
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where Co is the normahzation constant. Here, the eigenvalue f3 is an arbitrary complex 
number. 

It is easy to see that for any complex number f3 the state Eq.( p.2|) reduces to the 
coherent state \/3) = X]-^o "T^l'^) limit r/ — > 0, as expected. However, it 

does not have a number state limit for arbitrary f3 since number states are eigenstates 
of with non-negative integer eigenvalues. Further, we would like to have truncated 
states which are finite superpositions of the number states just as the binomial states 
are. With this in mind, we must choose (3 = ^/r]M, where M is a non-negative integer. 
In this case it is easy to see that the coefficients C„ are truncated 

{0, when n > M, 

(,r^y ^' ^Co, whenn<M. (^.3) 

Here the normalization constant CQ{ri,M) is obtained as 
Co(r/,M): 



Ln=0 



1 

2 \M 



Vj [(M-n)!]%!j Jm\Lm{-\^) 



(2.4) 



where X = J{1 — r])/-!] and Lm{x) is the Laguerre polynomial 



1 / M \ 

Lm{x) = Y.-[,, V-lTx-. (2.5) 



Inserting Eq. ( |2.3|) and Eq. (|2.4| ) into Eq. (|2.1| ), we obtain the desired solution ||?7, {(3 = 
^M)) = \\r^,M) 

1 Ml 
||r/, M) = , V A*^-'^ (2.6) 

which is a finite linear superposition of number states. 

We now consider the limiting cases of above state Eq. (|2.6|) as number and coherent 
states. First consider the limit 77 — 1. From the number-state expansion Eq. ( p.6|) , it 
follows that 

A*^-"M' 

^" = VMW.{M-n)\ ^'-^^ 
namely, Hr^, M) \M). Then, in the different limit — > 0, M — 5> 00 with ^JfjM = a a. 
real constant, we have 

^' M", A-"M"^a", Co exp(-aV2). (2.8) 



{M-n)\ 

and therefore Eq.( [^.6D reduces to the coherent state \a). 

The above discussion shows that the state \\ri, M) may be considered as an interme- 
diate state which interpolates between a number state and a coherent state. 
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2.2 Connection with photon-added coherent states 



The state Eq. (|2.6| ) can be written in more elegant form. By making use of \n) = -^lO), 
we can write Eq. (|2.6| ) as 



\V,M) 



M\Lm{-\^) 



J2 f^Vat)"A^"-" 



.n=0 



n 



|0) 



M!Lm(-A2; 



(2.9) 



where we have used the binomial formula. 

Furthermore, thanks to the equation (real A in our case) 

D{-X)a^D{X) = + A, 
where -D(A) is the displacement operator 

D(X) = exp A(a^ — a 
we can rewrite Eq. ( p.9|) in the form 

HM) = 



(2.10) 



(2.11) 



1 

m\lm{-x': 

D{-\)\X,M) 



D{-X)a^^ D{X)\0) 
D(-A)a^^|A) 



where |A) = L'(A)|0) is a coherent state and 

|A,M)e 



M\Lm(-X^) 



at^lA) 



(2.12) 



(2.13) 



is a so-called photon- added coherent state or excited coherent state 0. So from Eq. ( p. 131 ) 
we conclude that our new intermediate number-coherent states are displaced excited 
coherent states. 

However, we would like to point out that our states are very different from the 
photon-added coherent states. The photon-added states are an infinite superposition of 
number states from M to infinity, while our states are a finite superposition of number 
states from to M. 



3 Nonclassical Properties 

In this section we shall investigate the statistical and squeezing properties of ||?7, M) , 
with special emphasis on the comparison with those of the BS. 
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3.1 Photon statistics 



The easily- derived relation 



M{M - 1) ■ ■ ■ (M - A; + l)LM_fe(-A2) 



1/2 



\v,M-k), {3.1] 



ior k < M and a'^||?7, M) = for A; > M, gives the mean value of (A^) and (A^^) 

. M^M^, ,3.2) 

^^2^ ^ M(M - 1)Lm-2(-A^) + MLm-i(-A^) ^ ^3 3^ 

and Mandel Q-parameter [0 

M) ^ <M)! _ 1 . ,M - - M^-f^. (3.4) 

(iV) Lm-i(-A^) Lm(-A^) 

If Q{ri,M) < [resp. Q{ri,M) > 0), the field in the state \\ri,M) is sub-Poissonian 
(re5]).super-Poissonian). Q{ri,M) = corresponds to Poissonian statistics. 

For a fixed M, there are two extreme cases, r/ = (or A = cxd) and = 1 (or A = 0). 
It is easy to see that 

Qiv,M)^{;^' (3.5) 

in agreement with the Q-parameter values for number states and the vacuum state. Here 
we have used the fact -LAf(O) = 1 and Lm{x) / Ln{x) for m < n and x — oo. 

Fig. 1 is a plot of Q{r], M) with respect to rj for M = 2, 50, 100. The Q-parameter of 
the binomial states (BS) is also presented in the figure {Q = —i] for any M). From the 
figure we observe that the field in ||?7, M) is sub-Poissonian except for the case rj = 0. 

We say that a field is antibunched if the second-order correlation function ^''•^^(O) = 
(a) aa) / {a^ a)"^ < In fact, the occurrence of antibunching effects and sub-Poissonian 
statistics coincides for single mode, time-independent fields such as the state \\rj^M) of 
this paper. So the field H?], M) is antibunched except at the point r/ = 0. 

3.2 Squeezing properties 

Define two quadratures x (coordinate) and p (momentum) 

x = ^(a + a)^, p= (a - a'^) . (3.6) 
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Then we can easily obtain their variances (Ax)^ = (x^) — (x)^ and (Ap)^ = (p^) — (p)" 



(Ax)' = - + 



2 Lm(-A2) Lm(-A2) 



(3.7) 



where L^)(x) is the associated Laguerre polynomial defined by 

If (Ax)' < 1/2 (or (Ap)' < 1/2), we say the state is squeezedin the quadrature x (or p). 

Fig. 2 is a plot showing how the variance (Ax)' depends on the parameters t] and M. 
When ?7 = 0, (Ax)' = 1/2 since the state is just the vacuum state and in this case the 
field is not squeezed. Then, as t] increases the field becomes squeezed until maximum 
squeezing is reached; then the squeezing decreases until it disappears at a point r^o 
depending on M. We note that r/o < 1 when M > since (Ax') = M + 1/2 > 1/2 when 
— > 1. 

We also observe from Fig. 2 that the larger M, the stronger the squeezing, and the 
wider the squeezing range. 

It is known that the optimal signal-to-quantum noise ratio for an arbitrary quantum 
state 

(3.10) 

has the value 4Ns{Ns + 1) which is attainable for the usual coherent squeezed state0. 
For a coherent state the maximal ratio is 4A^s, where Ng is the mean value of the number 
operator for the quantum state. 

For the intermediate number-coherent state ||?7,M), the signal-to-quantum noise ra- 
tios for different parameters rj and M are shown in Fig. 3. The ratio for 77 = and rj = 1, 
which correspond to the vacuum state and number state respectively, is zero. For other 
1], we find from Fig. 3 (a) that the larger M, the larger the ratio. Fig.3 (b) gives plots of 
4(A^)((A^) + 1) ((A^) is given by Eq.(3.2)), 4(A^) and the ratio for the state \\r], M) with 
M = 10. We find that 

(1) the ratio for ||?7, M) is always smaller than the value 4{N){{N) + 1), which is in 
accord with the general result 0; 

(2) for some values of rj the ratio is larger than 4:{N). We observe that the states 
with ratio larger than 4(A^) correspond to squeezed states (see Fig. 2.). 
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4 Quasi-probability distributions 



Quasi-probability distributions [T^l in the coherent state basis turn out to be useful 
measures for studying the nonclassical features of radiation fields. In this section we 
study the Q-function (also called the Husimi function) and the Wigner function of the 

state \\ri, M). 

One can prove that (see Appendix A), if two states \ip)a and satisfy \il))a = 
D{a)\ilj), where D [a) = " " is the displacement operator, the Q and Wigner func- 
tions of \ip)a are simply a displacement of those of {ip), namely 

QiPh). = QiP - c^hh W^(/5)|^>. = WiP - a)|^). (4.1) 

So the Q-function and the Wigner function of the state ||?7, M) are easily obtained from 
those of the photon-added coherent states given in p 

= mv,M)\'= ■ (4.2) 

^.(,) . ?(zl)^^^L^±Mexp(-2|,P). (4.3) 

The Q-function Eq.( ^.2| ) has a 2M-fold zero at the position j3 = —A, which signals 
the nonclassical behaviouiQ. These zeros are related to the negative parts of the Wigner 
function, since the Q-function can be defined as a smoothed Wigner function. Fig. 4 
gives plots of the Wigner function of ||?7, M) for M = 3 and different t]. One can clearily 
see the negative parts, except for the case 1] = which corresponds to the vacuum state 
whose Wigner function is simply a Gaussian centered at the origin. As t] increases from 
0, the Gaussian distribution continuously deforms to the Wigner function of the number 
state 1 3). 

We can also study squeezing properties from the Q-function by examining the defor- 
mation of its contours. Fig. 5 is the contour plot of Q functions for M = 10 and different 
rj. We see that, when we increase rj, the contour is squeezed in the x direction until a 
maximum squeezing is reached. Then the contour deforms to the shape of a banana, 
which occupies a wider range in the x direction and the squeezing is reduced. Finally, 
we obtain a circular contour for larger t] corresponding to no squeezing (c.f. Fig. 2). 



We thank the referee for this remark. 
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5 Generation of intermediate states 



The main difference between the intermediate states described herein and photon-added 
coherent states is that the former are a finite superposition of number states. This 
suggests the possibility of an experiment to produce these states using the method 



proposed in [11 



We can also generate the state ||?7,M) by using the interaction of a photon and a 
two-level atom with an external classical driving field A in a cavity. In the rotating wave 
approximation, the Hamiltonian(/;, = 1) is 

H = Ho + V, 

I 

Ho = ujN + A{a^ + a) + -cJofTg, 

V = gi^a'a^ + aa^), (5.1) 

where = \e){e\ — \g){g\, cr+ = \e){g\ and cr_ = \g){e\ are atomic operators, g is the 
one-photon coupling constant, ujq and u are the atomic transition frequency and cavity 
resonant mode frequency respectively, and we take the driving field A to be real and 
constant. The interaction Hamiltonian is 

Hi{t) = Uo\t)VUo{t), Uo{t) = e-^-^«* = e-^/2'^*^-^^*('^^+»)e-^'"°*"^ (5.2) 

Using the following relation (see Appendix B) 

Uo\t)aUo{t) = e-''^'D{-A/uj)aD{A/io), (5.3) 

where D{A/uj) is the displacement operator, we have 

Hi{t) = gDi-A/cu) (e'("-"«)*aV_ + e-'("-"")*aa+) D{A/to). (5.4) 

Now we consider the on- resonance case, uj = ujq. Then the interaction Hamiltonian is 
t ime-indep endent 

Hi = gD{-A/uj){a^a_ + aa+)D{A/uj) (5.5) 
and therefore its time evolution operator is 

Uj{t) = e-'^'' = D{-A/u})e-'3'^''^''-+'"'+^D{A/uj). (5.6) 

Suppose that the field is initially prepared in the vacuum state |0) and the atom in the 
excited state |e); namely, at t = 0, the system is in the state |0) ® |e). At time t we have 

f//(t)|0) ® \e) = D(-A/cu)e-*f*("^"-+""+)D(A/cu)|0) ® |e). (5.7) 



When ^ft -C 1, we have 

Ui{t)\0)®\e) = \0) 0\e) -igt[D{-A/iu)a^D{A/uj)\O)]0\g). (5.8) 

If the atom is detected in the ground state \g), the field is reduced to the state ||?7, 1) 
with 7] = uo'^/{A^ + cj2). 

The state ||r7,M) (M > 1) can be generated by a multiphoton generahzation of the 
HamiltonianEq. ( |5. ![ ) , namely, V = g{a)^^(j^ + a'^^a^). 

Note that the parameter A depends on the external driving field and is a tunable 
parameter. In particular, for large enough M, we can control the output state to be 
either a number or a coherent state by tuning the parameter A. 

Finally we may infer the presence of these new intermediate states to first order in 
an idealized non-linear optics experiment. Consider a nonlinear Mach-Zehnder interfer- 
ometer with a Kerr medium in one arm. The output state is the displaced Kerr state 



D{i)UK{l)\\), Uk{i) ^ exp (^7«^V) , (5.9) 

where D{^) is the displacement operator and 7 = 2xL/v, L is the length of the Kerr 
medium, v the appropriate phase velocity inside the medium and x the third-order sus- 
ceptibility. When = —A, and 7 is small enough, the above states can be approximated 
as 

|0) + ^7A'||A,2) (5.10) 

showing the presence of the state \7], 2) in first order. In general, if we use a (25* + l)th- 
order nonlinear Kerr medium modelled in the interaction picture by [l^ 



HKerr = (5^^ = (^^j^^^ - 1) ■ ■ ■ {N - S) , (5.11) 

we can find ||?7, M) when 75 is small enough. 

6 Interaction with a two-level atomic system 

In this section we turn to the interaction of the state \\7],M) with a simple two-level 
system in the framework of the two-photon Jaynes-Cummings model (JCM) |jl^, |15 
described by the following Hamiltonian {h = 1) 

H = ua^a + -cuoda -|- g{a^'^(T^ + a^cr+) = Hq + V (6.1) 



with 



Hq = ua)a + -cuoCTs, V = g{a^^a^ + a^cr+). 
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The notation is as in Eq. ( ^.l]) , but now g is the two-photon couphng constant for 
transition \g) ^ |e). Suppose that, at the initial time t = 0, atom and field are decoupled 
and the atom is initially in the excited state |e), while the field is in the intermediate 
number-coherent state \\r]^M). Then the combined atom-field wave function at time t 



is obtained as 



M A A 

l^/W) = EC'n(^,M)[cos(5„t)-z— sin(5„t)]e*^*|e)® |n) 
n=0 25 

M Q 

-i -^Cniv, M) sm{6j)e-'^g) ® |n + 2), (6.2) 

n=0 



where 



f^n = (?\/(n + l)(n + 2), 6n = J^ + ni A = cuq - 2cu. (6.3) 



For simplicity, we only consider the on-resonance interaction case A = as in |T8 
whereupon Eq.(|6.2|) simphfies to 

M M 

\Mt)) = E Cniv, M) cos(fi„t)|e) ®\n)-tJ2 Cniv, M) sin(fi„t)|(7) ®\n + 2). (6.4) 

n=0 n=0 

We now discuss some quantum characteristics of the system arising from the equation 
Eq. 



6.1 Atomic population inversion 

Atomic population inversion is an important atomic observable in the JCM and is defined 
as the difference between the probabilities of finding the atom in the excited state and 
in the ground state. From Eq.( |6.4| ), the atomic population inversion is obtained as 

M 

W{t) = {as) = J2 |C„(r/,M)pcos(2fi„t). (6.5) 

n=0 

Fig. 6 gives the inversion vs. scaled time t = gt for different M and t]. From Fig. 6, we 
observe that the atomic population inversion exhibits the conventional Rabi oscillation 
for the M- number state limit [r] 1). In fact, in the limit rj ^ 1, Eq. ( |6.5| ) is simplified 
as 

W{t) = cos{2nMt) (6.6) 

with frequency 2fiM = 2^[(M+l)(M+2)]i/2 2MgioT high enough (A^), see Fig.6(a)). 
In the coherent state limit we observe the collapse-revival phenomenon, as we expect, 
with a revival time tcs which can be estimated as n/g ||T^ for high enough (A^) (that 
is, revival frequency Qcs = ^n/tcs ~ 2g) (Fig.6(d)). For the general intermediate case 
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(Fig.6(b, c)), remnants of both behaviour are seen; namely, an oscillation of frequency 
Qm modulated by the frequency Qcs with modulated amplitude dependent on the pa- 
rameter 7] and M. 

6.2 Field entropy 

We now consider the cavity field observables, beginning with entropy which is a measure 



of the amount of chaos or lack of information about a system [16|. The entropy S* of a 



quantum-mechanical system is defined as [17, 18 



^ = -Tr(pln(p)), 

where p is the density operator of the quantum system and the Boltzmann constant k is 
equal to unity. For a pure state, S = 0; otherwise S > 0, and it increases with increasing 
number of microstates with decreasing statistical weight. 

In this subsection we study the time evolution of the field entropy in our system. 



Barnett and Phoenix have proved that the field entropy Sj equals the atomic entropy 
Sa if the total initial state is a pure state. From Eq. (|6l^) the atomic reduced density 
operator pa can be easily obtained as 



where 



Pa = Trf{p) = Pii|^)(^| + Pi2|^)(e| +P2i|e)(^| +P22|e)(e|, (6.7) 



M 

Pii = ElCn(r?,M)|^in2(fi„t) 

n=0 
M 



P22 = ElCn(r/,M)|2cos2(n„t) (6.8) 

n=0 

M-2 

Pi2 = P21 = E Cn+2{v, M)Cn{r], M) cos(f2„+2t) sin(f2„t). 

n=0 

Then the field and atomic entropy Sa = — Tra(pa In(pa)) can be expressed as 

Sf = Sa = — vr+ ln(7r_|_) — 7r_ ln(7r_), (6.9) 

where 7i± are eigenvalues of the atomic reduced field density operator pa 

1 



7r± = -(1 ± ^(p22 - pii)2 + 4|pi2|2). (6.10) 

The field entropy Sf as a function of r is presented in Fig. 7. It is clear that Sf 
is a periodic function of time and it exhibits the conventional oscillation for the M- 
number state limit. As in the case of coherent initial states, the field entropy during 
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the time evolution is dynamically reduced to zero at revival time tji which means the 
cavity field can be periodically found in pure states, and reaches a maximum at and 
falls quickly to a minimum at r = 7r/4, 37r/4. Furthermore, for the general intermediate 
case, the field entropy has more minima as shown in Fig. 7 (b, c) due to the frequency 
modulation. 

6.3 Q-function 



The quasi-probabihty distribution Q-function is defined as |19 



1 



Qi/3) = -{P\pf\f3), 



TT 



where P/ = Tra(p) is the field reduced density operator 

M 

Pf = Cm{v,M)Cn{r],M)[cos{Qrnt)cos{Qnt)\n){m\ 

m,n=0 

+ sin(fi„t) sm{nj)\n + 2)(m + 2|] , 
and 1/3) is the coherent state. So the Q-function of the cavity field is 



(6.11) 



TT 



M 



Y^^Cn{7],M) COs(fi„t) 



n=0 



M 



E \ ^. Cniv, M) sin(l]„t) 



n=0 



n + 2)! 



(6.12) 

In Fig. 8 we give contour plots of the Q-function at different times r for r] = 0.1, 0.8. 
At time r = 0, the Q-function has only a single peak and the field is in the pure 
quantum state \ri, M) (c.f. Fig. 5). With the development of time, the Q-function begins 
to separate into two peaks. The smaller r/, the faster the seperation. At time r = 7r/2, 
the Q-function exhibits the most separation and the field is in a mixed state since the 
entropy reaches its maximum. Then two peaks begins to merge together and finally 
combine in a single peak at time r = vr, where the field is in a pure state with vanishing 
entropy. 

6.4 Photon number distribution 

The photon number distribution Pn{t) of the field described by the reduced density 
matrix p/ is given by 

P„(t) = {n\pf\n). (6.13) 
Inserting Eq. ( |6.11| ) into Eq. ( |6.13| ) we find the photon number distribution at time t 



Pn{t) = |C„(r7,M)|2cos2(^]„t) + |C„„2(r7,M)|2sin2(fi,_2t). 



(6.14) 
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Fig. 9 shows the behaviour of the photon number distribution at times r = 0,7r/4, 
7r/2, 37r/4 and vr. From these figures we can observe that the photon number distribution 
exhibits strong oscillation at time r = 7r/4 and 3/47r for the intermediate states. In fact, 
at those times, the field is a superposition of two components (see Fig. 8) and its emtropy 
decreases rapidly to a minimum (see Fig. 7). Partial interference between two component 
results in strong oscillation of the photon number distribution. However, the oscillation 
is not perfect (see below). Nevertheless it is perfect at the slightly earlier time r = 7r/4— ^ 
(see dashed lines in Fig.9(a,b)). 

This effect is not hard to understand. In fact, at r = 7r/4, we have the following 
approximate result (see Apendix C) 



Pn{t) 



|C„(r/,M)|2sin^ 



1 

n I r 

2 



(6.15) 



r/2(M-n + 2)2(M-n + 1)2_ 
for the high enough (A^) case. Eq.( |6.15D is a strongly oscillating function which explains 



the large oscillations of the photon number distribution. However, due to the additional 
term r/2 = tt/8, the function sin^[(n — l/2)r]^=7r/4 cannot be zero for any integer n; in 
other words, the oscillation is not perfect. However, Pn{t) is zero at the slightly earlier 
time r = it/A — ^, where ^ is chosen to make (n — l/2)r a multiple of vr. 

From Fig. 9(c) we also observe that the photon number distribution at r = vr is simply 
a displacement by 2 from that at the time r = 0. For the large photon number case, 
this fact can be proved analytically. Using Eq. ( |C.4| ) in Appendix C, we have 

sin(^]„_2^) ~ sin(n7r - n/2) = (-l)"+\ 
cos{Qnt) ~ cos[(n + l)n + n/2] = 0. 
So the photon number distribution Eq. (|6.14|) at r = tt becomes 

Pn{7f/g) = |a-2(^,M)|2 = P„_2(0). (6.16) 
In the same way we find that, at r = 7r/2, the photon number distribution is 

Pn{n/2g) = ^ (|C„(r/,M)|2 + |C„_2(r7,M)|2) = i (p„(0) + P„(7r/(7)) , (6.17) 

namely, the average of the photon number distributions at r = and r = vr. In Fig. 9(c) 
this fact can be clearly observed. 

7 Conclusion 

In this paper we have described new states \\ri,M) which interpolate between number 
and coherent states and have investigated their various properties. Unlike photon-added 
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coherent states, to which they are related, these states are a finite superposition of num- 
ber states. We also analyzed the limiting cases rj ^ 1 and rj ^ 0, M ^ oo corresponding 
to number and coherent states respectively. Salient statistical properties of \\ri, M) such 
as the sub-Poissonian distribution, the anti-bunching effect and the squeezing effects 
were investigated for a wide range of parameters. The non-classical features of thse 
states for certain parameter ranges were demonstrated in terms of the quasiprobability 
distributions, the Q and Wigner functions. We also proposed an experiment to generate 
these states, inferring their presence in certain non-linear systems. 

We then considered the interaction of these interpolating number-coherent states with 
a two-level atomic system, exemplified by the two-photon Jaynes-Cummings Model. We 
first studied the dynamics of atomic population inversion. On an intuitive level, one 
expects that the response of the atomic system will vary between the Rabi oscillation 
typical of an initial number state, and the collapse-revival mode for an initial coherent 
state; and indeed this is what one obtains. We found that it exhibited the conventional 
Rabi oscillation for the M-number state limit with frequency VLm (~ '^Mg) and the 
collapse-revival phenomenon for the coherent state limit with rivival frequency Qcs ~ 251. 
For the general intermediate case, remnants of both behavior were seen; namely, an 
oscillation of frequency modulated by the frequency Q^^s with modulated amplitude 
dependent on the parameter 77 and M. 

We further investigated the field observables, the entropy, Q-function and photon 
number distribution. It is of interest that the photon number distribution exhibit strong 
oscillation at r = 7r/4, 37r/4. At those times, the field entropy falls rapidly to a minimum 
and the Q-function separates into two peaks, which means that the field is a superpo- 
sition of almost pure states and interference between components of the superposition 
state leads to strong oscillation of the photon number distribution. An approximate 
analytical solution is presented to explain this result. 

The remarkable properties of these intermediate number-coherent states provide a 
useful tool for theoretical investigation of model systems; their generation by non-linear 
systems tempts us to believe that the states found in this paper may play an important 
role in quantum optics. 
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Appendix A: Displaced quasiprobability distributions 

For the Q-function, we prove Eq. ( [4.1D as follows: 

Q(/3)l^)^ = mD{am\' = mD{-P)Diam\' 

= \{p - a\^)\' = QiP - a)\^}, (A.l) 

where we have used the relation 

D{5)D{'y) = D{6 + 7)e^('^^*-T"^*) = D{6 + 'yy^'^^^^'\ (A.2) 

for arbitrary complex numbers S and 7. From the following definition of the Wigner 
function [pO| 



r) 00 

W{P) = -Y.{P,k\p\p,k), (A.3) 

where \P,k) = D{i3)\k) = e^"'''~^*°'\k) is the displaced number state {\k) is the number 
state) and p = \ri, M){r], M\ is the density matrix of the states considered, we can prove 
the second relation in Eq. (|4.1|) in the same way as in the Q-function case. 

Appendix B: Proof of formula Eq . ( [5t3| ) 



In this appendix we give a proof of Eq. ( ^.31 ). We use the following formula 

00 ^_^y„ 
n=0 



-^Ge^ = Y: [F, [F, • ■ ■ , [F, G] ■■■]]. (B.l) 



n copies 

For the case in hand 



F = -iutN -itA{a^ + a), G = a. (B.2) 
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It is easy to see that 



[F, G] = iuta + lAt, 
[F,[F,G]]=iut[F,Gl 

[F, [F, [F, G]]] = zut[F, [F, G]] = {zutYlF, G], 



^ ^ 

n copies 

= M)"a + M)"A/cj = {iujt)''D{~A/uj)aD{A/uj), (B.3) 



where D{A/uj) is the displaced operator. Substituting Eq. ( p.3|) into Eq. ( p.l| ) we obtain 
the formula Eq.( |5.3|) . 



Appendix C: Photon number distribution for large 
photon number 

In this appendix we present an analytical treatment of the photon number distribution 
in the large photon number regime. The photon number distribution of the two-photon 
JC model with initial state |e) C?> J2n Cn\n) can be obtained as 



Pn{t) = |C„,r cos" i^^in + l){n + 2)rj + |C„_2r sin' (^V - ' (^-1) 

where r = gt is the scaled time as before. 

Here we only consider an initial field state which is narrower than that of a coherent 
state. For a distribution {|C„p} we can calculate the variance as 

(n - nf = {N^) - {N)\ (C.2) 

For the coherent state we have (n — n)^ = n. So for highly excited coherent states 
where n —>■ cxd, we have n n. In the following we only consider a distribution {|C,ip} 
narrower than the Poisson distribution, namely 

{n-nf<n. (C.3) 

So for large enough n we also have n n. In this case, we have 



3 3/1 



'n + l)(n + 2) ^ ^/ri^ + 3n = n\ll + - = n + -= [n- - ] + 2 



n — l)n ^/-n? — n = nJ 1 ~ ^ ~ 2' (^•^) 



Futhermore, when r = 7r/4 we have 



cos 













j = cos^ 






= sin^ 


(n ^ r 






' 2. 




A 2/ . 



. (C.5) 



Substituting Eq. (|C.5| ) into Eq. (|C.l| ) we obtain the approximate photon number distri- 
bution at r = 7r/4 

r 



Pn(t) = (|C„p + |C„_2p) sin^ 



n r 

2 



(C.6) 



From Eg. ( |C.6| ) we find that, for the initial field whose photon distribution is narrower 
than a Poisson distribution, the photon number distribution at r = 7r/4 exhibits strong 
oscillation. However, P„(7r/4) cannot be zero for any n due to the term r/2 = tt/S and 
the oscillation is not perfect. Nevertheless, the oscillation is perfect at a slightly earlier 
time T = tt/A — ^, as indicated in Ref. |[T8[| (for initial coherent state) and Fig.9(a,b). 

For coherent states we further have |C„_2p ~ jC'^ip. So the photon number distri- 
bution is 



P„(t) =2e-"— sin^ 
n! 



n 



(C.7) 



which is just the result given in Ref. [jT8|. 

Now we turn to the analytical approximate result Eq.( |6.15| ). For the intermediate 
state, we can write the variance as 



{n - umY = {N)m-i{N)m - {N)l + {N)m, 



(C.8) 



where hm = {N)m = (f], M\N\t], M). In general, we have (A^)a/-i < {N)m < M. 
For large enough um, or M, we have {N)m-i ~ {N)m and therefore {n — umY ~ ""-m 
which leads to n ~ um- So the result Eg. ( |C.6|) is valid for the intermediate state case. 
Furthermore, the distribution |C„(r7, M)p and |C„_2(?7, M)p are related by 

(1 — 7])'^n{n — 1) 



\Cn-2{r],M)Y 



(C.9) 



r72(M-n + 2)2(M-n + 1)2' 

Sunstituting Eq.( p.9| ) into Eq.( |C.6| ), we finally obtain Eq.( |6.15D . 

Eq. ( |C.4| ) can also be used to explain the behaviour of the photon number distribution 
at r = 7r/2 and vr (see Eqs.( CT , |6T7D and Fig. 9). 
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Figure 1: Mandel's Q parameter for M=2, 50, 100. 





Figure 3: The signal-to-quantum noise ratio for Hr], M): (a) The ratio for different M; 
(b) Comparison of p , 4(iV)((iV) -f- 1) and 4(iV) for M = 10. 
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Figure 4: Wigner function of ||77, M) for M 




3 and rj — 0.1, 0.4, 0.7 and 1. a — x + iy. 
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Figure 7: Entropy of the field as a function of scaled time r. (a) M = 4 and r] — 0.9999 
(the initial field state is the number state |4)); (b) M = 70 and 77 = 0.8; (c) M = 70 and 
T] = 0.1; (d) M = 200 and rj = 0.005. 



24 




25 




35 40 45 50 „ 55 60 65 70 



0- 



(c) 


Ti=o.8 nn 


Thick solid line: x=0 




Thin solid Une: t;=7c/2 




Dashed line: x=t 




Tl=0.1 

In \^ 







35 40 45 50 „ 55 60 65 70 



Figure 9: Number distribution of the photon field at different times r = 0,7r/4,7r/2,37r/4 
and TT for 77 = 0.1, 0.8. In (a) and (b) we also present the distribution at slightly eariher 
time T — ^, where ^ is chosen as 1/140 and 1/180 for 77 = 0.1 and 77 = 0.8 respectively. 
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